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Thermoelectric  devices  can  provide  clean  energy  conversion  and  are  environmentally  friendly.  However, 
few  studies  have  been  conducted  to  understand  the  dynamic  characteristics  of  (thermoelectric  gener¬ 
ators)  TEGs.  The  stability  of  output  performance  plays  a  very  important  role  in  TEG.  Based  on  the 
thermoelectric  effect  and  heat  transfer  theory,  a  dynamic  model  for  waste  heat  recovery  in  a  general  TEG 
was  developed  to  asses  the  influence  of  heat  reservoir  and  heat  sink.  The  validity  of  the  theoretical  model 
was  demonstrated  by  experiments  and  it  was  found  that  the  model  prediction  agrees  well  with  the 
experimental  results.  The  dynamic  response  characteristics,  such  as  hot  and  cold  semiconductor  surface 
temperatures,  maximum  output  power  and  system  efficiency,  were  studied  using  this  dynamic  model. 
The  results  showed  that  enhancing  heat  dissipation  on  cold  side  is  the  crucial  to  improve  the  output 
performance  of  TEG,  and  that  the  fluctuation  of  hot  reservoir  leads  to  rapid  change  of  output  power 
generated  by  TEG,  which  is  dangerous  to  electric  devices.  The  dynamic  model  proposed  in  the  present 
study  can  be  used  in  the  design  and  operation  analysis  of  TEGs. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

As  energy  shortage  and  environmental  deterioration  are  grow¬ 
ing,  energy  saving  and  emission  reduction  have  become  global 
obligation  and  responsibility.  Besides  being  environmentally 
friendly,  the  thermoelectric  technology  has  many  other  advantages 
such  as  being  highly  reliable,  adapting  to  different  kinds  of  heat 
reservoirs  and  having  no  moving  parts.  Therefore,  the  generation 
technology  has  received  more  and  more  attention  in  scientific 
community,  especially  in  the  aspect  of  waste  heat  recovery  [1]. 

In  order  to  improve  the  energy  conversion  efficiency  of  TEG,  re¬ 
searchers  have  completed  a  number  of  experimental  and  modeling 
studies.  Niu  et  al.  [2]  setup  a  low-cost,  simple  configuration  TEG  unit 
with  commercially  available  Bi2Te3  based  on  the  (thermoelectric 
module)  TEM  for  an  anticipated  maximum  power  generation  of 
about  150-200  W.  Casano  and  Piva  [3]  analyzed  TEG  on  the  basis  of 
the  experimental  data  for  ‘open’  and  ‘closed’  circuit  voltage,  electric 
power  output  and  conversion  efficiency  as  a  function  of  temperature. 

In  the  case  of  simulated  study,  a  numerical  model  with  taking 
contact  effects  and  actual  structure  into  account  was  developed  by 
Liang  et  al.  [4].  It  was  found  that  the  matched  load  of  TEG  was  not 
equal  to  the  internal  resistance.  Hsiao  et  al.  [5]  applied  a  mathe¬ 
matic  model  of  thermoelectric  module  with  applications  on  waste 
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heat  recovery  from  automobile  engine  to  establish  fundamental  an¬ 
alyses  of  TEM.  Zhou  [6]  presented  a  two-stage  TEG,  analyzed  and 
optimized  the  performance  in  the  objective  function  of  efficiency  and 
output  power.  And  Chen  et  al.  [7]  analyzed  the  effects  of  design  factors 
on  the  performance  of  the  two-stage  TEG,  and  they  presented  a  novel 
configuration  for  a  two-stage  thermoelectric  heat  pump  system 
driven  by  a  two-stage  TEG  [8].  The  effect  of  the  heat  transfer  irrever¬ 
sibility  on  the  performance  of  TEG  was  also  researched  [9,10].  In 
consideration  of  temperature  dependence  of  thermoelectric  prop¬ 
erties  and  external  heat  transfer  irreversibility,  Chen  et  al.  [11] 
established  a  new  model  of  multi-element  thermoelectric  generator. 
The  results  show  that  the  temperature  dependence  of  thermoelectric 
properties  has  a  significant  effect  on  the  power,  efficiency  and  optimal 
variables  of  TEG.  A  linear  phenomenological  heat  transfer  model  [12] 
and  a  radiative  heat  transfer  model  [13]  were  also  presented  by  Chen 
et  al.  Hsu  et  al.  [14]  deliberated  the  issue  how  to  transform  the  energy 
of  waste  heat  into  electricity  for  vehicles  and  motorcycles. 

In  addition,  some  authors  not  only  have  conducted  experimental 
studies,  but  also  developed  simulating  models.  Crane  and  Jackson 
[15]  focused  on  simultaneously  optimizing  both  the  heat  exchanger 
geometry  and  the  thermoelectric  geometry  for  optimal  net  waste 
heat  recovery  performance.  In  our  previous  study,  we  have  also 
conducted  an  optimization  on  low-temperature  waste  heat  ther¬ 
moelectric  generator  system  by  experiment  and  modeling  [16]. 

Most  of  the  researchers  focused  on  the  static  properties  of  TEGs, 
and  few  studies  on  the  dynamic  characteristics  have  been  carried 
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Nomenclature 

Subscripts 

a 

air 

A 

area  (m2) 

b 

base 

cp 

specific  heat  capacity  (J  kg-1  K-1) 

a/b 

material 

E 

electromotive  force  (V) 

c 

cold,  channel 

h 

convective  heat  transfer  coefficient  (W  m-2  K_1) 

ce 

ceramic  substrate 

H 

height  (m) 

ch 

characteristic  length 

I 

electric  current  (A) 

con 

contact 

I< 

thermal  conductance  (W/K) 

conv 

convective  heat  transfer 

l 

length  (m) 

cu 

copper  conducting  strip 

n 

number 

ef 

effective 

P 

perimeter  (m),  power  (W) 

ex 

heat  exchanger 

Q 

heat  flux  (W/m2) 

f 

fin 

Qv 

inner  heat  generation  per  unit  volume  (W/m3) 

fc 

flow  channel 

a 

heat  transfer  rate  (W) 

h 

hot 

Qv 

inner  heat  generation  (W) 

i 

inner 

r 

electric  resistance  (Q) 

L 

load 

R 

thermal  resistance  (K/W) 

m 

matched 

S 

surface 

max 

maximum 

T 

temperature  (I<) 

N 

N-type  semiconductor 

u 

velocity  (m/s) 

P 

P-type  semiconductor 

W 

flow  rate  (m3/h) 

sc 

semiconductor 

y 

coordinate 

si 

solder  layer 

ss 

stainless  steel 

Greek  letters 

w 

water 

a 

Seebeck  coefficient  (V/I<) 

d 

thickness  (m) 

Abbreviation 

V 

efficiency  (%) 

Nu 

Nusselt  number 

X 

thermal  conductivity  (W  m-1  K-1) 

Pr 

Prandtl  number 

P 

density  (kg/m3) 

Re 

Reynolds  number 

V 

Kinematic  viscosity  coefficient  (m2/s) 

th 

hyperbolic  tangent  function 

T 

time  (s) 

A 

difference 

out.  Moreover,  the  TEG  output  performance  is  susceptible  to  heat 
reservoir.  The  running  process  is  an  extremely  unstable  process  in 
practice,  when  taking  waste  heat  as  the  hot  reservoir  of  a  TEG; 
meanwhile,  the  stability  of  output  performance  plays  a  very  impor¬ 
tant  role  in  TEG.  Therefore,  it  is  essential  to  study  the  TEG  dynamic 
response  characteristics  along  with  changes  of  heat  reservoirs. 

In  this  study,  a  theoretical  dynamic  model  of  a  thermoelectric 
generator  for  waste  heat  recovery  has  been  built.  At  the  same  time, 
a  corresponding  thermoelectric  device  is  constructed  to  validate 
the  model.  Based  on  the  dynamic  model,  the  TEG  dynamic  response 
characteristics  have  been  modeled  when  temperatures  and  flow 
rates  of  heat  reservoirs  follow  step  changes. 

2.  System  modeling 

According  to  Seebeck  effect,  the  TEG  open-circuit  voltage  de¬ 
pends  mainly  on  the  temperature  difference  between  the  cold  and 
hot  sides  of  semiconductor.  That  is  to  say,  changes  of  the  temper¬ 
ature  difference  lead  to  changes  of  the  open-circuit  voltage,  and  the 
TEG  output  performance  will  change  as  well.  Therefore,  the 
unsteady-state  heat  transfer  model  is  built  to  solve  the  TEG  tem¬ 
perature  distribution  first.  And  the  TEG  dynamic  model  is  setup 
after  the  corresponding  performance  calculation  model  is  built. 

2.1.  Unsteady -state  heat  transfer  model  [17] 

To  build  the  unsteady-state  heat  transfer  model,  it  is  essential  to 
conduct  an  unsteady-state  the  heat  transfer  analysis  first.  Similar 
with  Ref.  [17],  the  basic  configuration  for  the  TEG  dynamic  study 


follows  the  general  design  illustrated  in  Fig.  1.  Thermoelectric 
modules  are  sandwiched  between  heat  exchanger  and  flow  chan¬ 
nel  with  hot  fluid  that  can  be  treated  as  a  convective  heat  transfer 
system. 

A  TEM  consists  of  two  ceramic  substrates,  a  certain  number  of 
semiconductors,  corresponding  solder  layers,  copper  conducting 
strips  and  thermal  insulation  material,  as  illustrated  in  Fig.  2.  The 
semiconductors  are  connected  thermally  in  parallel,  but  electrically 
in  series. 

For  a  well-designed  thermal  insulation  packing  TEM,  the  heat 
leakage  through  the  TEM  surrounding  could  be  neglected. 
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1.  ceramic  substrate  2.  semiconductor  3. solder  layer 
4.  copper  conducting  strip  5.  thermal  insulation  material 


Fig.  2.  A  thermoelectric  module. 


Furthermore,  both  the  specific  heat  capacity  and  the  volume  of  air 
in  TEM  are  so  small  that  the  heat  leakage  through  the  surrounding 
of  TEM  parts  could  also  be  neglected.  In  addition,  the  temperature 
differences  between  TEM  parts  are  small  because  temperature  of 
hot  reservoir  studied  is  not  high.  So  the  radiation  heat  transfer  in 
TEM  could  be  neglected.  Therefore,  the  TEG  model  can  be  processed 
in  a  one-dimensional  heat  transfer  way,  as  shown  in  Fig.  3. 

According  to  heat  conduction  theory,  the  one-dimensional  un¬ 
steady  heat  conduction  equation  on  the  hypothesis  of  physical 
parameter  being  constant  can  be  shown  as  [18]: 

a  t  ,a2T 

pCpfr  ~  Xdyi  +  qv  (1) 


The  corresponding  boundary  conditions  are  provided  to  solve 
the  equation  (2).  As  shown  in  Fig.  3,  the  boundary  conditions  on  the 
cold  and  hot  sides  of  heat  exchanger  belong  to  the  third  boundary 
condition.  The  boundary  conditions  of  other  center  sides  in  TEG 
which  are  regarded  as  the  interrelated  and  complementary  con¬ 
ditions  belong  to  the  second  boundary  condition.  There  are  the  left 
and  right  boundary  faces  on  each  center  side,  including  Sce,h,  5cu,h, 

Ssl,h>  SSCih,  $sc,o  Ssl.c*  Scu.c  and  Sce, o  which  are  expressed  as  S+  h,  S~e  h, 
c+  c-  c+  c-  c+  c-  c+  c-  c+  c-  c+  ’  c-’ 

°cu,h’  °cu,h’  Jsl,h»  °sl,h»  Jsc,h’  Jsc,h’  Jsc,c»  Jsc,c»  Jsl,c’  Jsl,c»  Jcu,c*  Jcu,c* 

S+  cand  S~e  c  respectively.  According  to  the  treatment  of  Peltier  heat 
and  Joule  heat,  the  boundary  condition  of  each  part  in  TEG  is  listed 
as  following. 

(1)  On  hot  side  of  flow  channel  (SfC,h) 

—  fttAr,fc  (^h  —  ^fc,h)  (3) 

x=sfc,h  v  7 

(2)  On  hot  side  of  ceramic  substrate  (Sce,h) 


,  .  ar 

-Afc/lfc— 


ax 


,  ar 

-/exnex/iext)— 


x=St, 


,  A 

—  —^ceElce^ce-r- 

ax 


x=S~ 


(4) 


where  p,  cp  and  A  are  density,  specific  heat  capacity  and  thermal 
conductivity  respectively,  t  and  y  are  time  and  axis  coordinate 
respectively,  and  qv  is  inner  heat  generation  per  unit  volume. 

Unlike  common  power  generation  technology,  the  thermo¬ 
electric  generator  involves  three  basic  effects,  Seebeck  effect,  Pelt¬ 
ier  effect  and  Thomson  effect.  Besides,  there  are  two  accessorial 
effects,  Joule  effect  and  Fourier  effect.  Seebeck  effect  generates 
electromotive  force.  Peltier  effect,  Thomson  effect  and  Joule  effect 
lead  to  Peltier  heat,  Thomson  heat  and  Joule  heat,  respectively. 
Therefore,  the  TEG  inner  heat  generation  includes  Peltier  heat, 
Thomson  heat  and  Joule  heat,  which  are  conducted  an  analyzed  in 
the  following. 

Firstly,  Peltier  heat  generates  on  the  two  sides  of  semiconductor, 
so  it  can  be  disposed  as  the  boundary  heat  flux.  Secondly,  Thomson 
heat  is  small  enough  to  be  neglected  since  it  is  a  second-round 
effect.  Finally,  Joule  heat  can  be  regarded  as  flowing  equally  to 
the  two  sides  of  conductors  [19],  so  it  can  also  be  disposed  as  the 
boundary  heat  flux.  In  conclusion,  the  TEG  inner  heat  generation  is 
equal  to  zero.  The  TEG  unsteady-state  heat  transfer  equation  (1)  is 
translated  into 


ar 

pCpfo 


(2) 


(3)  On  hot  side  of  copper  conducting  strip  (Scu, h) 


(5) 


(4)  On  hot  side  of  solder  layer  (5si,h) 


(5)  On  hot  side  of  semiconductor  (5sc,h) 


,  .  ar 

-AslHslAsi— 


1  t2 

ynslrsi  = 


.  dT 

— Ascns  cAsc^ 

+  nScttITs,h 


x=S~ 


1  j2 

—  —I  riscTsc 


(7) 


skJ»smM.sm>.sgA.sKM  -  hot  side  of  flow  channel,  ceramic  substrate,  copper  conducting  strip,  solder  layer,  semicondictor 
-  cold  side  of  heat  exchanger,  ceramic  substrate,  copper  conducting  strip,  solder  layer,  semiconductor 


Fig.  3.  One-dimensional  heat  transfer  model  of  thermoelectric  generator. 
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(6)  On  cold  side  of  semiconductor  (SSCtC) 


,  A  dT 
—AscTIscAsct^ 


1 


dT 


x=Sfc 


—I  Hsc^sc  +  tlsc^/rsc,c  —  — ^sI^sI^sIq^ 

1  j2 

-Jns\r si 


X=S7C 


(8) 


(7)  On  cold  side  of  solder  layer  (SsitC) 


,  .  dT 

*~As\ns\As\fc 


1 


dT 


x=St 


^  ns\rs\  —  — ^CuHcu^CUq^ 


x=S7i 


1  r2 

—  2  tlcu^CU 

(9) 


(8)  On  cold  side  of  copper  conducting  strip  (Scu,c) 


in 


(16) 


where  Pf  and  Af  are  perimeter  and  area  of  cross  section  of  fin 
respectively,  and  hc  is  coefficient  of  convective  heat  transfer  on  cold 
side  of  heat  exchanger.  While  hh  in  equation  (3)  is  coefficient  of 
convective  heat  transfer  on  hot  side  of  flow  channel.  The  two  co¬ 
efficients  are  solved  by 


(17) 


where  /Ch  is  characteristic  length,  and  Nu  is  Nusselt  number. 
Because  the  structures  of  the  convective  heat  transfer  are  different, 
the  equation  of  Nusselt  number  of  hot  side  of  flow  channel  is  dif¬ 
ferent  from  that  of  cold  side  of  heat  exchanger. 

For  internal  pipe  flow  of  hot  fluid,  the  Nusselt  number  is 
expressed  best  by  Gnielinski  correlation  [20] 


,  ar 

-XcutlcuTlcuTr- 

dX 


X—ScUC 


1  t2  i  a  dT 

+  2^  ^cufcu  =  —  ^cetlce^ce^ 


(10) 


X=Sr n 


(9)  On  cold  side  of  ceramic  substrate  (Sce,c 


,  dT 

-AceftceAzeQ^ 


,  dT 

—  —  ^extlexTipxT- 


x=S+iC 


ax 


x=S- c 


(11) 


Nl,  C/7 8)(Reh  -  1000)Prh 
1  +  12.7v478(Pr2/3  -  1  j 

Reh  =  2300- 106 


(18) 


where  Reh  and  Prh  are  Reynolds  number  and  Prandtl  number  of  hot 
reservoir  respectively,  kc  and  /Ch,  fc  are  length  and  characteristic 
length  of  flow  channel  respectively,  and  /  is  Darcy  resistance  coef¬ 
ficient  of  turbulent  tube  flow.  ZCh  fC,  and /are  given  respectively 
by 


(10)On  cold  side  of  flow  channel  (5fCtC) 


,  A  dT 

^ex71ef  ex  ^ 


X — Se> 


—  hcAe f  ex  (Tex, c  —  Tc) 


(12) 


where  AfC,  Ace,  Acu,  hb  he  and  Aex  are  thermal  conductivities  of  flow 
channel,  ceramic  substrate,  copper  conducting  strip,  solder  layer, 
semiconductor  and  heat  exchanger,  respectively.  Ace,  Acu,  As\,  Asc 
and  Aex  are  cross  sectional  areas  of  ceramic  substrate,  copper  con¬ 
ducting  strip,  solder  layer,  semiconductor  and  heat  exchanger, 
respectively;  nce,  ncu,  ns i,  nsc  and  nex  are  numbers  of  ceramic  sub¬ 
strate,  copper  conducting  strip,  solder  layer,  semiconductor  and 
heat  exchanger,  respectively;  rcu,  rs\  and  rsc  are  electric  resistances 
of  copper  conducting  strip,  solder  layer  and  semiconductor, 
respectively;  Th  and  Tc  are  temperatures  of  hot  reservoir  and  cold 
reservoir,  respectively.  TfC,h  is  hot  side  temperature  of  flow  channel, 
lex,  c  is  cold  side  temperature  of  heat  exchanger,  and  Aefjc  and  Aef>ex 
are  the  effective  areas  of  flow  channel  and  heat  exchanger  which 
are  expressed  as  respectively 


Jch,fc 


Pfc 


(19) 


Reh 


*Vch,fc 


(20) 


/  =  (1.821gReh-1.64)-2  (21) 

where  PfC  and  AfC  are  perimeter  and  area  of  cross  section  of  flow 
channel  respectively,  Vh  is  kinematic  viscosity  coefficient  of  hot 
reservoir,  and  Uh  is  velocity  of  hot  reservoir  which  is  expressed  as 


(22) 


where  Wh  is  the  flow  rate  of  hot  reservoir. 

For  the  external  cooling  air  flow,  a  fined  heat  exchanger  is 
showed  in  Fig.  4.  The  expression  of  the  Nusselt  number  can  be 
calculated  by  Ref.  [21] 


^ef,fc  —  tlceTlce 


(13) 


^ef,ex  —  ^ex,b  +  VfAexf 


(14) 


where  Aex,b  and  Aex,f  are  surface  areas  of  heat  exchanger  base  and 
fin  respectively,  and  rjf  is  fin  efficiency  of  heat  exchanger.  And  the 
fin  efficiency  is  given  as 


m 


th(mHf) 

niHf 


(15) 


where  th  is  hyperbolic  tangent  function,  Hf  is  the  height  of  fin,  and 
m  is  a  constant  which  is  calculated  by 


Nuc  =  0.664 Rel^Pr^3,  Rec  <  5  x  105(Laminar  flow)  ^3) 
Nuc  =  0.037 Ret/5Prl/3,  Rec>  5  x  105(Turbulent  flow) 

where  Rec  and  Prc  are  Reynolds  number  and  Prandtl  number  of  cold 
reservoir,  respectively.  Rec  is  given  by 


Rec 


wc/ch,ex 

Vc 


(24) 


where  Zc h,  fc  is  the  characteristic  length  of  heat  exchanger  which  can 
be  taken  as  the  height  of  fins  [17],  vc  is  kinematic  viscosity  coeffi¬ 
cient  of  cold  reservoir,  and  uc  is  velocity  of  cold  reservoir,  which  is 
expressed  as 
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heat  source  Th 


heat 

leakage 

KX(T^- 


„  KciT^~T<) 

- ►  cold  source  Tc 

Fig.  4.  Structure  diagram  of  heat  transfer  system  of  thermoelectric  generator. 


hot  side  of  semi 
-conductor  Tscc 

f 

'  rj2  f 

\  oclTxb 

r~)  c§> 


1 1 

'  rJ2  \ 

\ 

|  alTxc 

cold  side  of  semi 
-conductor  Tscc 

Uc 


wc 

Af 


(25) 


where  Wc  is  flow  rate  of  cold  reservoir,  and  Af  is  the  passageway 
area  between  fins. 


2.2.  Performance  calculation  model 


The  TEG  output  performance  is  evaluated  by  calculating  output 
power  and  the  system  efficiency  [2].  According  to  Seebeck  effect, 
the  electromotive  force  of  thermoelectric  circuit  is  expressed  by 

E( t)  «  na[TSC5h(T)-rsc,c(T)]  (26) 

where  Tsc>h  and  Tsc>c  are  the  hot  and  cold  side  temperatures  of 
semiconductor  respectively,  n  is  the  number  of  thermoelectric 
couple,  and  a  is  Seebeck  coefficient.  The  current  of  thermoelectric 
circuit  is 


/(t) 


Eft) 

E[  +  rCon  +  fL 


(27) 


Qsc,h  —  (^h  ^sc,h) 


(29) 


Qsc, C  =  n/Cc(TSc,c  -  Tc) 


(30) 


where  J<b  is  the  thermal  conductance  from  hot  reservoir  to  hot  side 
of  semiconductor,  and  I<c  is  the  thermal  conductance  from  cold 
reservoir  to  cold  side  of  semiconductor.  K h  and  J<c  are  given  by 


I< h 


1 

\/Kh  +  VKfc  +  VI<ce  +  1/Kcu  +  1/Ksi  +  Rc  on 


(31) 


Kc 


1 

1/KC  +  1/Kex  +  1/^ce  +  Wen  +  Wsl  +  ^con 


(32) 


where  K h  is  thermal  conductance  from  hot  reservoir  to  hot  side  of 
flow  channel,  I<c  is  thermal  conductance  from  cold  reservoir  to  cold 
side  of  heat  exchanger,  JCfC,  I<ce ,  I<cu,  I<s\  and  Kex  are  thermal  con¬ 
ductance  of  flow  channel,  ceramic  substrate,  copper  conducting 
strip,  solder  layer  and  heat  exchanger  respectively,  and  Rcon  is  the 
TEM  contact  thermal  resistance. 

In  addition,  Peltier  effect  shows  up  as  heat  absorption  on  hot 
side  and  heat  release  on  cold  side  of  semiconductor,  and  Joule  effect 
shows  up  as  heat  release  on  hot  and  cold  sides  of  semiconductor.  As 
showed  in  Fig.  4,  Qsc.h  and  Qsc.c  are  also  given  by 


Qsc,h  —  ft 


(33) 


Qsc,c  —  ft 


Mlsc,c  +  ^  vsc  '  rsc,s 


V  DC, 11 


(34) 


where  rsc  and  I<sc  are  electric  resistance  and  thermal  conductance  of 
a  couple  of  semiconductors  respectively,  rsiiSc  is  contact  electric 
resistance  between  solder  layer  and  semiconductor  on  the  hot  side, 
and  rsc,si  is  contact  electric  resistance  between  semiconductor  and 
solder  layer  on  the  cold  side. 

The  heat  loss  is  neglected  when  the  TEG  heat  transfer  model  is 
taken  as  one-dimensional  model.  Based  on  energy  conversation 
law,  the  heat  semiconductor  releases  to  cold  reservoir  and  the 
output  power  add  up  to  the  heat  semiconductor  absorbs  from  hot 
reservoir: 


Qsc,  h  —  Qsc,  c  + 

Combine  (25)  to  (35), 


(35) 


where  n,  rcon,  and  rL  are  inner  electric  resistance,  contact  electric 
resistance  and  load  electric  resistance,  respectively.  So  the  output 
power  is  given  by 


Ei-ft)  =  l\ 


n2a2[Tsch(T)-TSCtC(z)]2 

- ; - ,2 - rL 

(n  +  Icon  +  rL) 


(28) 


Arsc 


T  t  —  (jh  ~  Tc)(n  +  rcon  +  ijj 

sc,h  sc.c  -  Bn2a2  +A(r.  +  rcon  + 


where 

A  =  1  +  (jH+i-c)Ksc 


(36) 


(37) 


When  load  electric  resistance  is  equal  to  the  matched  electric 
resistance,  the  output  power  reaches  maximum.  So  the  matched 
electric  resistance  is  solved  to  obtain  the  maximal  output  power. 
The  expression  of  the  matched  electric  resistance  is  deduced  in  the 
steady  condition  as  follows. 

Fig.  4  presents  semiconductor  absorbing  heat  from  hot  reservoir 
and  releasing  heat  to  cold  reservoir.  According  to  thermodynamic 
theory,  semiconductor  can  regarded  as  heat  engine  working  be¬ 
tween  hot  reservoir  and  cold  reservoir,  which  is  the  dotted  part  in 
Fig.  4.  Based  on  Fourier  Law,  the  heat  which  semiconductor  absorbs 
from  hot  reservoir  and  releases  to  cold  reservoir  is  expressed  as 


n_Th  , 

I<C  Kh 

Substitute  (36)  to  (28) 

n2a2rL(Th  -  Tc)2 


Pi  = 


[Bn2  a2  +  (n  +  rcon  +  rL)Ay 


(38) 


(39) 


The  equation  (39)  is  taken  the  derivative  with  respect  to  load 
resistance.  Then  let  the  derivative  equation  be  equal  to  zero, 
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^  =  o 
arL 

Therefore,  the  matched  load  resistance  can  be  obtained 
Bn2a 2 

rm  =  H  +  rCOn  H  * — 


Substitute  (41)  to  (28),  the  maximal  output  power  is 

n2«2(>\2 n  +  ABn2«2)  [TSC;h(T)  -  Tsc,c(t)]2 


<(t)  = 


(2Arj  +  Bn2  a2) 


(40) 


(41) 


(42) 


the  air  cooling  is  used  on  the  cold  side.  To  enhance  the  heat  dis¬ 
sipation,  six  heat  exchangers  are  fixed  on  the  cold  side  of  the 
modules.  At  the  same  time,  two  fans  are  installed  on  the  heat 
exchanger  for  the  forced-convection  heat  transfer. 

In  addition,  there  are  eight  thermocouples  used  in  data 
acquisition  system.  As  shown  in  Fig.  5,  two  thermocouples  are 
located  at  the  inlet  and  outlet  of  flow  channel,  and  other  thermo¬ 
couples  are  located  averagely  on  both  surfaces  of  thermoelectric 
module  in  one  side  of  flow  channel.  In  order  to  measure  the  tem¬ 
perature  in  real  time,  a  series  of  virtual  instruments  from  National 
Instruments  Corporation  (NI)  are  employed  here. 


According  to  Newton’s  law  of  cooling,  the  TEG  heat  absorption 
capacity  is 

Qh(T)  =  hhAfc[Th(x)-TeXth(z)}  (43) 

Finally,  the  TEG  system  efficiency  is  calculated  by 

(T)  Pmax(T)  n2«2 (A2r,-  +Afin2«2)  [Tsc,h(t)  -  rsc,c(t)]2 

QhM  Me  f  (2 M  +  B"2«2)2  ph(T)  -  TeXih(  T)] 

The  performance  calculation  model  is  built  by  deducing  the 
expressions  (42)  and  (44).  So  the  TEG  dynamic  model  has  been 
setup  with  the  addition  of  the  unsteady-state  heat  transfer  model. 

3.  Experimental  validation 

3.1.  Experimental  setup 

An  experimental  setup  which  is  designed  to  validate  the  TEG 
dynamic  model  has  been  established.  The  physical  and  schematic 
diagrams  of  the  setup  are  shown  in  Figs.  5  and  6,  respectively.  The 
system  consists  of  a  thermoelectric  converter,  a  hot  fluid  loop, 
a  cooling  device  and  a  data  acquisition  system. 

In  the  hot  fluid  loop,  a  single  stainless  steel  channel 
(600  mm  x  100  mm  x  18  mm,  and  2  mm  wall  thickness),  which  hot 
water  flows  through,  is  lined  on  each  side  with  eighteen  thermo¬ 
electric  modules  (TEC1-12708-NET1,  a  kind  of  thermoelectric 
cooling  module  with  127  thermoelectric  couples,  and  40  mm 
length,  40  mm  width  and  3.5  mm  height).  All  the  modules  are 
connected  electrically  in  series.  Hot  water  that  is  used  to  simulate 
waste  heat  is  supplied  by  the  circulating  constant  temperature  bath 
(temperature  range:  room  temperature  -  100  °C).  A  water  circu¬ 
lating  pump  is  located  to  monitor  the  flow  rate  of  hot  water.  And 


Fig.  5.  Physical  diagram  of  the  experimental  setup  of  TEG  applied  in  waste  heat 
recovery. 


3.2.  Comparison  of  simulation  and  experimental  results 

Experiments  are  conducted  under  the  specific  operating  con¬ 
ditions  as  follows:  the  temperature  of  hot  water  is  between  50  °C 
and  90  °C,  the  air  temperature  is  constant  28.5  °C,  and  the  flow 
rates  of  hot  and  cold  fluid  are  2.76  m3/h  and  43.2  m3/h,  respectively. 
Besides,  the  physical  parameters  and  dimension  parameters  of  each 
part  in  TEG  are  listed  in  Tables  1  and  2. 

Out  of  flatness  of  the  contact  surface  between  TEM  and  the 
convective  heat  transfer  system  that  is  caused  by  the  installment 
of  thermocouples,  the  contact  thermal  resistance  is  so  large  that  it 
has  to  be  considered  into  the  model.  The  contact  thermal  resist¬ 
ance  can  be  gained  by  experiment,  and  be  replaced  by  the  cor¬ 
responding  thermal  resistance  of  grease  with  approximately 
8.6  x  10~4  m  of  thickness  and  1.5  W  m-2  K-1  of  thermal  con¬ 
ductivity  in  the  simulating  process.  Meanwhile,  the  contact  elec¬ 
tric  resistance  needs  to  be  calculated  for  the  reason  that  the 
number  of  thermoelectric  modules  is  large  enough,  even  if  a  single 
module  has  a  small  contact  electric  resistance.  The  matched  load 
resistance  can  be  gained  by  experiment.  So  the  whole  TEG  contact 
electric  resistance  can  be  calculated  by  the  expression  (41 ).  Finally, 
the  calculated  value  of  the  contact  electric  resistance  is  approx¬ 
imately  6  Q. 

To  compare  simulated  results  with  experimental  data,  the 
simulation  must  have  the  same  operating  conditions  as  the 
experiment  besides  based  on  an  accuracy  system  model.  Each 
experiment  is  repeated  at  least  three  times  to  obtain  reproducible 
results.  And  the  mean  vale  of  the  three  effective  experimental  re¬ 
sults  acts  as  the  final  experimental  data.  Here,  the  TEM  ceramic 
surface  temperatures  and  output  power  are  regarded  as  the  study 
objects.  These  comparisons  are  shown  in  Figs.  7  and  8. 

Fig.  7  shows  the  temperature  variations  on  the  hot  and  cold 
ceramic  surfaces  with  different  hot  fluid  temperatures  in  steady 
states.  The  steady  states  mean  the  TEG  temperature  distribution 
remains  constant  with  time.  Furthermore,  because  Peltier  effect 
and  Joule  effect  reach  the  maximum  in  a  short  circuit,  they  have  the 
maximum  impact  on  TEG  temperature  distribution  right  now.  So 
Fig.  7  is  gained  in  a  short  circuit.  As  it  can  be  seen  in  Fig.  7,  the 
experimental  results  nearly  coincide  with  the  simulation  results. 
The  experimental  temperatures  on  the  hot  side  are  a  little  lower 
than  the  simulated  temperatures.  And  the  maximum  value  of  var¬ 
iation  is  within  2  °C. 

Output  power  versus  load  resistance  at  three  different  hot  fluid 
temperatures  is  plotted  in  Fig.  8.  All  the  simulated  results  are  a  little 
higher  than  the  experimental  data.  However,  the  relative  error 
between  experimental  and  simulated  results  is  within  20%.  At  hot 
fluid  temperature  of  90  °C,  the  error  is  only  1.9%. 

Therefore,  the  steady-state  simulated  results  and  experimental 
data  shown  in  Figs.  7  and  8  indicate  not  only  a  very  similar  trend, 
but  also  acceptable  relative  errors.  It  is  demonstrated  that  the 
simulated  results  are  in  good  agreement  with  the  experimental 
data  in  the  steady-state  operation,  so  the  reasonability  of  TEG  dy¬ 
namic  model  is  partially  verified. 
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I - Circulating  constant  temperature  bath  2-Circulating  pump  3-Flow  channel  4-Thermal  insulating  layer 
5 -Thermoelectric  module  6-Heatsink  7-Axial  fan  8-K-type  thermocouple  9-Ammeter  10-Voltmeter 

I I-  Load  resistance  box  12-Switching  power  supply  13-Data  acquisition  system  14-Desktop  computer 


Fig.  6.  Schematic  diagram  of  the  experimental  setup  of  TEG  applied  in  waste  heat  recovery. 


4.  Dynamic  simulation 

Based  on  the  established  dynamic  model,  the  dynamic  response 
characteristics  of  surface  temperatures  of  semiconductor,  max¬ 
imum  output  power  and  system  efficiency  are  studied  when  tem¬ 
peratures  or  flow  rates  of  heat  reservoir  generate  step  changes.  The 
initial  condition  in  which  TEG  operates  can  be  proposed  as  follows: 
hot  water  temperature  is  60  °C,  air  temperature  is  30  °C,  and  flow 
rates  of  water  and  air  are  1.38  m3  h-1  and  21.6  m3  h-1,  respectively. 

4.1.  Dynamic  characteristic  with  step  change  of  heat  reservoir 
temperature 

Fig.  9  presents  the  dynamic  response  characteristics  of  hot  and 
cold  semiconductor  surface  temperatures  with  step  changes  of  heat 
reservoir  temperatures.  There  are  four  step  changes  in  total  as 
shown  in  Fig.  9.  To  analyze  the  influence  of  each  step  change  on  TEG 
clearly,  the  set  value  of  the  time  interval  between  two  adjacent  step 
changes  must  assure  each  temperature  in  TEG  reaches  at  least  99% 
of  the  steady-state  value.  And  the  TEG  state  is  called  as  the  quasi¬ 
equilibrium  state  right  now.  Here,  it  is  found  that  the  time  inter¬ 
val  can  be  set  as  400  s  to  reach  this  state  by  simulating.  The  first 
step  change  is  hot  water  temperature  jumping  from  60  °C  to  70  °C 
at  the  time  of  zero.  Secondly,  hot  water  temperature  follows  a  step 
change  from  70  °C  to  60  °C.  Thirdly,  air  temperature  follows  a  step 
change  from  30  °C  to  20  °C.  Lastly,  air  temperature  returns  to  30  °C 
in  the  way  of  step  change. 

As  shown  in  Fig.  9,  the  hot  semiconductor  surface  temperature 
increases  nearly  to  the  steady-state  value  in  a  surprisingly  short 


Table  1 

Basic  physical  parameters  of  each  component  in  thermoelectric  generator. 


Parts 

Specific  heat 
capacity  J/(  kg  K) 

Thermal 
conductivity  W/ 
(m-K) 

Density 

kg/m3 

Electrical 
resistivity  O  m 

Semiconductor  158 

1.5 

7300 

1  x  10"5 

Solder  layer 

167 

33 

8560 

5  x  lO"7 

Copper 

386 

398 

8930 

2  x  10"8 

conducting 

strip 

Ceramic 

765 

36 

3975 

- 

substrate 

Flow  channel 

644 

12 

7600 

- 

Heat 

871 

162 

2660 

- 

exchanger 

time,  and  the  cold  semiconductor  surface  temperature  increases 
slowly  compared  with  the  hot  semiconductor  surface  temperature 
at  the  first  stage.  At  the  second  stage,  the  variations  are  identical 
but  the  change  directions  are  opposite  to  those  at  the  first  stage  for 
the  hot  and  cold  semiconductor  surface  temperatures.  Both  the  hot 
and  cold  semiconductor  side  temperatures  display  an  exponential 
decay,  but  the  decrease  of  the  cold  side  temperature  is  bigger  than 
that  of  the  hot  side  temperature  at  the  third  stage.  The  relationship 
of  the  response  characteristics  between  the  fourth  and  the  third 
stage  is  just  like  that  between  the  second  and  the  first  stage. 

Since  the  thermal  resistance  from  hot  reservoir  to  hot  side  of 
semiconductor  is  very  small,  hot  side  temperature  of  semi¬ 
conductor  possesses  the  characteristic  of  increasing  rapidly  at  the 
first  stage.  The  main  reason  may  be  due  to  that  the  thermal 
resistance  of  semiconductor  occupies  a  considerable  proportion  in 
the  whole  thermal  resistance.  Based  on  the  trends  of  the  hot  and 
cold  semiconductor  surface  temperatures,  the  trend  of  the  tem¬ 
perature  difference  between  the  two  surfaces  shown  in  Fig.  9  is 
then  easy  to  understand.  In  addition,  the  response  time  of  the  hot 
semiconductor  surface  temperature  at  the  first  stage  is  much 
shorter  than  that  of  the  cold  semiconductor  surface  temperature  at 
the  third  stage.  According  to  heat  transfer  theory,  this  phenomenon 
can  be  explained  by  the  much  larger  the  convective  heat  transfer 
coefficient  in  hot  reservoir  than  that  in  cold  reservoir. 

Fig.  10  shows  the  response  characteristics  of  maximum  output 
power  and  system  efficiency  at  the  same  operating  conditions  with 
Fig.  9.  It  can  be  seen  that  maximum  output  power  reaches  a  peak  in 
a  short  time,  then  decreases  slowly  to  a  stable  value  at  the  first 
stage.  It  is  clearly  shown  in  Figs.  9  and  10  the  whole  tendency  of 


Table  2 

Dimensions  of  each  component  in  thermoelectric  generator. 


Parts 

Cross  section 
length  mm 

Cross  section 
width  mm 

Height 

mm 

Number 

Semiconductor 

1.47 

1.47 

1.3 

4  x  2  x  127 

Solder  layer 

1.47 

1.47 

.05 

4  x  2  x  127 

Copper 

4 

1.6 

.4 

4  x  127 

conducting 

strip 

Ceramic  substrate  40 

40 

.65 

4 

Flow  channel 

100 

18 

600 

1 

Heat  exchanger 

80 

80 

6 

1 

base 

Heat  exchanger 

80 

4 

38 

8 

fin 


’  In  table  means  the  parameter  needn’t  to  be  used  in  calculation. 


208 


X.  Gou  et  al.  /  Energy  52  (2013)  201-209 


Time  (s) 


Fig.  7.  Experimental  data  and  simulating  results  for  the  hot  and  cold  ceramic  surface  Fig.  9.  Hot  and  cold  semiconductor  surface  temperatures  versus  time  with  step 
temperatures  versus  hot  water  temperature  in  the  steady-state  operation.  changes  of  heat  reservoir  temperatures. 


maximum  output  power  are  consistent  with  that  of  the  tempera¬ 
ture  difference  in  hot  side  and  cold  side  of  semiconductor.  This 
result  can  be  explained  by  the  expression  (42). 

Though  the  heat  input  of  the  thermoelectric  converter  will 
change  when  the  heat  reservoir  temperatures  make  step  changes, 
its  relative  variation  is  much  lower  than  that  of  output  power. 
Hence,  system  efficiency  variation  depends  mainly  on  the  variation 
of  output  power.  The  reason  causing  the  oscillations  is  that  the  heat 
transfer  rate  is  finite  so  that  hot  side  temperature  of  flow  channel 
cannot  respond  as  quickly  as  hot  water  temperature. 

4.2.  Dynamic  characteristic  with  step  change  of  heat  reservoir  flow 
rate 

The  dynamic  response  characteristics  of  semiconductor  surface 
temperatures,  maximum  output  power  and  system  efficiency  with 
step  changes  of  heat  reservoir  flow  rates  are  described  in  Figs.  11 
and  12.  Each  step  change  of  heat  reservoir  flow  rates  takes  place 
at  the  same  moment  as  the  corresponding  step  change  of  heat 
reservoir  temperatures.  The  first  two  change  steps  are  that  flow 
rate  of  hot  water  makes  a  round  trip  between  1.38  m3  h_1  and 
2.76  m3  h-1.  In  this  process,  flow  rate  of  air  follows  the  third  step 


Fig.  8.  Experimental  data  and  simulating  results  for  output  power  versus  load 
resistance  for  three  different  hot  water  temperatures  in  the  steady-state  operation. 


change  from  21.6  m3  h-1  to  43.2  m3  h-1,  which  returns  to 
21.6  m3  fr1  at  the  fourth  step  change.  It  is  found  that  the  curves  in 
Figs.  11  and  12  have  the  same  trend  as  those  in  Figs.  9  and  10, 
respectively.  Accordingly,  adding  flow  rate  of  hot  water  could 
generate  a  similar  impact  with  increasing  its  temperature  on  TEG; 
while  for  cooling  air,  adding  flow  rate  has  a  similar  impact  on 
decreasing  its  temperature. 

It  can  be  seen  in  Fig.  12,  though  both  flow  rates  of  water  and  air 
have  changed,  the  variation  of  maximum  output  power  and  system 
efficiency  at  the  last  two  stages  are  larger  than  those  at  the  first  two 
stages.  This  result  shows  enhancing  heat  dissipation  on  cold  side 
generates  a  greater  improvement  on  TEG  output  performance. 

Finally,  combining  the  above  figures,  it  can  be  seen  that  the 
whole  tendencies  of  maximum  output  power  and  system  efficiency 
are  mostly  consistent  with  that  of  the  corresponding  temperature 
difference  between  hot  side  and  cold  side  of  the  semiconductor.  On 
this  view,  TEG  output  performance  is  mainly  depended  on  the 
temperature  difference  between  hot  side  and  cold  side  of  the 
semiconductor.  When  the  temperature  difference  increases,  TEG 
output  performance  can  be  improved.  In  addition,  when  hot  water 
temperature  or  flow  rate  change  rapidly,  output  power  generates 
a  sharp  change  remarkably.  Therefore,  in  order  to  protect  electric 
equipment,  it  needs  to  avoid  or  reduce  the  variation  of  output 
power  generated  by  fluctuation  of  hot  reservoir  in  system  design 
and  operation. 


Fig.  10.  Maximum  output  power  and  system  efficiency  versus  time  with  step  changes 
of  heat  reservoir  temperatures. 
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b.  In  order  to  protect  electric  equipment,  it  is  necessary  to  avoid 
or  reduce  the  variation  of  output  power  generated  by  fluctua¬ 
tion  of  hot  reservoir  in  system  design  and  operation. 

Therefore,  these  dynamic  simulation  results,  gained  by 
employing  the  dynamic  model  to  analyze  different  operating  con¬ 
ditions,  could  provide  guidelines  for  performance  optimization  and 
safe  operation  of  thermoelectric  generator.  The  results  have  pro¬ 
vided  a  useful  guideline  for  designing  TEGs. 
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Fig.  11.  Hot  and  cold  semiconductor  surface  temperatures  versus  time  with  step 

changes  of  flow  rates.  References 


E 

<u 


Fig.  12.  Maximum  output  power  and  system  efficiency  versus  time  with  step  changes 
of  flow  rates. 


5.  Conclusion 

In  this  study,  a  theoretic  dynamic  model  of  a  general  thermo¬ 
electric  generator  with  finned  heat  exchangers  is  presented,  and 
a  waste  heat  recovery  thermoelectric  generator  experimental  setup 
has  been  constructed  to  validate  this  model.  The  good  agreement 
between  the  steady-state  simulated  results  and  experimental  data 
has  approved  the  reasonability  of  the  system  dynamic  model. 
Therefore,  this  system  model  can  be  used  in  dynamic  simulation. 

By  means  of  step  changes  of  heat  reservoir  temperatures  and 
flow  rates,  the  dynamic  response  characteristics  of  temperatures  on 
the  hot  and  cold  semiconductor  surfaces,  maximum  output  power 
and  system  efficiency  are  studied.  The  TEG  output  performance  has 
different  dynamic  characteristics  at  different  operating  conditions. 
There  are  two  conclusions  obtained  as  follows: 

a.  Compared  with  enhancing  heat  transfer  on  hot  side,  enhancing 
heat  dissipation  on  cold  side  generates  a  greater  improvement 
on  TEG  output  performance. 
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